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Abstract 

The electric and magnetic self-fields can exist in the rotating plasma. 
The self-sustained rotation can be established in the plasma. The disturbed 
distribution function of rotating plasma is derived from the Vlasov equation. 
The propagation of waves in rotating plasma differs from that in the usual 
plasma. New terms of the Landau damping appear. The local rotational 
behaviour may become prevailing. 

1 INTRODUCTION 

Phenomena dealt with are not too rare. However, the causes of the phenomena 
(forces, fields) are rather weak and can often be supposed to be negligible. 

It is well known that many astronomical and laboratory objects represent 
rotating plasma. Therefore, the study of rotating plasma is of great theoretical 
and practical interest. 

It is assumed that plasma is neutral on the average. However, electric and 
magnetic fields can exist in plasma, because there are forces which oppose the 
electric one^'^. 

It is well known that in a wide variety of plasma devices plasma travels 
across the magnetic field at rates much higher than those explained by collisional 
processes. It has been assumed that the electron transport across the magnetic 
field arises from the E x B drift {E x B instability^). 

The radial electric field exists to balance out the v x B and centrifugal force 
terms^. 

The plane sheet stability of a surface charge has been considered for externally 
applied and mutually perpendicular fields {E _L B) which produce an equilibrium 
particle flow in the sheet^. The calculation was extended to a cylindrical shell in 
the inhomogeneous magnetic field^ (the size of a shell is identical to the particle 
gyro-orbit). 

The equilibrium properties induce a self-electric field for a charge slab or a 
cylinder (the fluid approximation'''^). 

The effect of equilibrium radial electric field on the trapped - particle stability 
has been investigated ^. The stabilization can result from the configuration 
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in which the radial potential assumes a hollow profile (the equilibrium radial 
electric field points inward). The theory also predicts the possibility of additional 
instability, when the E x B rotation is present (the rotationally driven trapped 
- particle mode). 

The development of a radial electric field Ep due to toroidal rotation has 
been studied for the toroidal geometry The ion distribution function becomes 
Maxwellian with temperature uniform on the magnetic surface. The poloidal flow 
decays. The density is distributed over magnetic flux surfaces by the Boltzmann 
factor with the effective potential energy that is the sum of a centrifugal potential 
and an electrostatic potential required for charge neutrality. The usual small 
gyroradius expansion of the Fokker-Planck equation was carried out, and it was 
shown, that the effects of centrifugal and Coriolis forces on the particle motion 
affect the transport. 

The rotational equilibrium in cylindrical, non-neutral plasma with an axial 
magnetic field has been analysed in detail 

The poloidal rotation of plasma due to the E x B drift has been studied in 
tokamaks (the so-called H — mode of improved confinement) with a radial elec- 
tric field and strong axial magnetic field. The E x B drift rotates the plasma 
poloidally (when the radial electric field exists). In non-i^-modc plasmas the 
poloidal rotation frequency is typically much less than the toroidal rotation fre- 
quency in tokamaks. The toroidal geometry of tokamaks becomes important in 
some theories of poloidal spin-up^^. The poloidal rotation instability (resulting 
from the poloidally asymmetric diffusion) leads to poloidal velocity shears which 
may quell the microturbulence. The nonlinear interplay between the poloidal 
spin-up and turbulence-driven anomalous transport leads to bifurcated equilib- 
ria (the L to H transition mode in tokamaks). The tokamaks can spontaneously 
develop a poloidal velocity shear (the "P-S" flows). The spin-up persists even 
under conditions where the growth rate is sonic 

The rotation is well known to induce new poles in the perturbed distribution 
function. The so-called resonant diochotron instability ^^'^^ involves the Landau 
damping at the resonant frequency uj = IQ (I is an integer, ft is the rotation 
frequency) . The role of Landau damping in cross-field electron beams and in an 
inviscid shear flow has been considered for the cylindrical geometry (perturba- 
tions with kz = 0) with unperturbed angular velocity u;o(r) = ^no,0 < r < b 
and oJo{r) = ^no{b/r)'^ ,b < r < a (no is a constant). A formal analysis of the 
resonance has been given in the Laplace transforms context The case of a 
perfectly conducting wall (at the impedance Z{ijj) = 0) was studied, and a slight 
shift in the eigenfrequencies from their locations for an induced Z{uj) ^ was 
found. The Landau damping can overcome the resistive destabilization (if there 
exist gradients in the electron density) for a strictly decreasing electron density 
in a cylindrical system, so that only one resistively unstable normal mode can 
exist - the fundamental one. 

The diochotron instability of a thin (^^^zzli ^1), tenuous {up <^uJBe)i cylin- 
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drical layer (a single-component plasma) has been investigated in a strong uni- 
form constant magnetic field using the Vlasov equation The particle gyro- 
radius was of the order of the mean radius of a layer. The theoretical model does 
resemble the geometry of some experiments which rely on the high energy parti- 
cles injection. Two distribution functions were considered, which give the same 
density {p{r) ~ [6(r — ri) — 0(r — r2)], where 6 is the Heaviside step function), 
but are different in the velocity space (the velocity spread is zero, or the parti- 
cles oscillate around the mean radius). The self-electric field of the equilibrium 
density (which is strictly radial) was included in the treatment (but not the self- 
magnetic field). The stability against electrostatic perturbations with frequency 
around —luose (where / is the asimuthal wave number of perturbations) and 
the infinite wavelength in the z-direction has been investigated. The perturbed 
distribution function was found, and the macroscopic sheared flow was shown to 
be not a relevant criterion for instability by itself. 

The gyrokinetic integral equations, used for studying the ion temperature 
gradient (ITG)- driven mode in a toroidal geometry at low plasma pressure, 
have been extended to include the equilibrium ion sheared flows The parallel 
and perpendicular sheared flows are destabilizing and stabilizing mechanisms, 
respectively. The gyrokinetic integral code was used to explore the sheared flows' 
effects on the ITG mode in a sheared slab geometry The mode structure and 
eigenfrequencies predicted by the integral code differ from the results derived by 
the differential approach (for higher radial eigenmodes). 

The further development of the above mentioned works and some new ideas 
are presented in this work. The purposes of the paper are: 

1) to show the separation of charges and to obtain rigorous expressions for electric 
and magnetic self-fields in rotating plasma; 

2) to demonstrate the self-sustained rotation origin; 

3) to derive the perturbed distribution function and to find the location of poles 
(the terms in the Landau damping); 

4) to propose the local-rotational description of plasma. 

The consideration is self-consistent in its character (from the first principles) 
and leans upon the microscopic properties of a matter, rather than upon the con- 
tinuous medium (fluid) models. The description can be applied to free rotating 
plasma objects (astronomical, for example) with some arbitrary shape and com- 
position, to plasma devices with cylindrical symmetry in particular (the arbitrary 
geometry is included as the limits of integration only) . 

The peculiarities, which differentiate this work from the cited papers, are as 
follows: 

- Not only the radial electric self-field Ep{p, z), but the magnetic self-field Bz{p, z) 
and smaller self-fields Ez{p-, z) and Bp(p, z) are taken into consideration (note that 
all these components are functions not only of p but of z as well). This consid- 
eration is necessary, since the behaviour of particles changes qualitatively: the 
solution depends not only on the small ratio | , but on arbitrary ratio of character- 
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istic velocities as well. All quantities (density, charge density, local composition, 
etc.) are not externally given but can be evaluated in a self-consistent manner. 

- The possibility of a steady-state rotation of plasma is shown. The self-sustained 
rotation of initially neutral plasma origins in the external magnetic field Bq as a 
consequence of the E x B drift {E is the self-field). 

- The perturbed distribution function is found (rigorously) proceeding from the 
Vlasov equation. The approach differs in carrying through the extra terms due 
to A:^, a self-consistent radial centrifugal force, the self- fields and in allowing the 
wave amplitudes to be space-dependent. 

- The new poles are shown to exist, and the location of these poles is found (the 
new terms in the Landau damping appear). Since the wave propagation in ro- 
tating plasma (an accelerated medium) differs from that in ordinary plasma, the 
local-rotational description of plasma is supposed to be possible (in the strong 
turbulence case, for example). 

2 THE ELECTRIC AND MAGNETIC FIELDS IN 
ROTATING PLASMA 

Rotating plasma possesses some peculiarities as compared to rotating objects 
which consist of neutral particles. We start with the rigid rotation for a simple 
explanation only (the substitutions for an arbitrary rotation are given below). 
The existence of the centrifugal force 

F = mOV, (1) 

(where (} is the angular frequency of the system, m is the particle mass, p is 
its distance from the axis of rotation) causes different effects on particles of 
different masses. As a result of these different effects, the p-dependence of particle 
concentration is bound to be different for particles of different masses. 
According to the Boltzmann distribution, 

na = noa exp j ' (2) 

where T is the system temperature, k is the Boltzmann constant, is the 
particle concentration on the axis of rotation, [/„ is the potential energy of a- 
sort particles. In the case of neutral particles 

(/„ = (3) 

The na{p) dependences for hydrogen plasma are shown schematically in Fig.l. 
Here i?o is the size of a system (the distance from the axis of rotation); symbol 
a is either H (hydrogen) or e (electron). 
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However, plasma consists of charged particles. Since the particles of differ- 
ent masses in plasma have different charges, this partial separation of particles 
produces a partial separation of charges. The negatively charged region is bound 
to lie near the axis of rotation, whereas the positively charged region lies near 
the system boundary (Rq). The electric field Eo{p,z) (in the polar coordinate 
system: the Q, direction is the z direction) exists as a result of charge separation. 
This field opposes the considerable separation of charges. 

As illustrated in Fig.l, there is some distance from the axis of rotation R'i{z), 
where the local charge equals zero. For < p < R'i{z) the plasma is negatively 
charged on the average; for R'i{z) < p < Rq the plasma is positively charged. 

The partial separation of charges gives rise to the p-dependence of charge 
density; that is, the plasma possessing the given charge density moves round a 
circle of definite radius. Therefore, circulating currents are inside the rotating 
plasma despite the fact, that all particles revolve with Cl on the average has a 
distinct direction) . For < p < R'l the current is opposite to the plasma rotation 
direction (since the charge is negative), whereas for R[ < p < Rq the current 
flows in the direction of rotation (since the charge is positive). In the general case 
the magnetic actions of currents are not compensated, and the magnetic field B 
exists. This is the magnetic field of a solenoid, whose axis coincides with the axis 
of system rotation. 

For rotating, fully ionized plasma (near macroscopic equilibrium) , which con- 
sists of elements with atomic number Ni, it follows in polar coordinates p, (p, z {z 
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axis coincides with the ft direction), that ne,ni can be taken from (2) with 

Ue = TT^ + e / Eop{p', zi)dp' + e / Eo,{pi,z')dz'+ 

^ Jo Jo 

r ^np'dp' - r ^^Pidz' , (4) 
Jo c Jo c 

MO? /-pi rzi 

Ui = 2^ - eNi Eop{p', zi)dp' - eNi Eo.ipi, z')dz'- 

/ ^Qp'dp' + / ^npidz' , (5) 
Jo c Jo c 

where nc{p, z) is the electron concentration, ni{p, z) is the ion concentration, — e 
is the electron charge (e > 0), m is the electron mass, M is the ion mass, Tg = 
Ti = T are electron and ion temperatures, r is the radius- vector: r = pep + ze^. 
The quantities riQa can be found from conditions: jfy^ Uadv' = NaKa, where 
Ka is the total number of a-sort particles, A^'g = 1 and K^, = NiKi + A, A is the 
electron surplus (for non-neutral plasma A 7^ 0). 

The Eo{p,z) components can be obtained from the system of integral equa- 
tions: 

Eop{p, z) = e [ '^^^l^f^\ i[p - picos{ipi - ip)]dpidzid(pi , (6) 

J(V) itT 



where 



q{pi,zi) = noiexp i-^Uij - noeexp i--^Ue\, (7) 



Ri = \J{z- ziY + p'^ + pj- 2ppicos{(fi - if) 

Eoz{p,z) = e[ ^(Pi' ^1) _ zi)pidpidzid<fi , (8) 
J{V) ill 

V is the system volume (its shape is arbitrary). 

Radius R'l can be found from the equation ne(Ri) = ni{R'i) using (2), (4)- (8). 
It follows from the Biot-Savart law, that the magnetic field is 

c J(y) R\ 

where c is the speed of light, dvi is the region with currents, Ri is the distance 
of this region to the point of observation, ji is the current density. 

The magnetic field due to rotation can be obtained from the following ex- 
pressions: 

1 /■ _ ^,)p,dp,dzidipi , (9) 

C J(v) R^ 



- / ^^\ i[picos{(fi -ip) - p]dpidzid(pi , (10) 

c J(v) Ri 
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where 

Mpi, zi) = enpiq{pi,zi), (11) 

with substitutions q{pi,zi) from (7), Uq, from (4), (5) and Eq{p',z') from (6), 
(8); the magnetic field near the axis of rotation at the system boundary is Bz 
from (10) iov z = L, p = 0, where 2L is the system size along the axis of rotation. 

The following substitutions need to be done for rotating plasma of complex 
composition: 

— ^ Y^ij^'ip , riQa — ^ nQ^; in the case of unsteady rotation 

n^p^ fp 

' Jo ^{p'^z)pdp- 

For astronomical objects the effects can be considerable; in this case the 
gravitational force plays important role. Note, that in the general case the field 
E^pe^p exists for free rotating objects with some arbitrary (non-symmetric) shape. 
In this case all quantities depend on coordinate as well. The additional terms 

rip rip 

e / EoipPidcpi , - eNi / Eoippidtpi 
Jo Jo 

appear in (4) and (5), respectively. The component Eq^ is taken from 

Eoip{p, z,ip) = e / pism{ip - ipi)dpidzidipi. 

J{V) -rti 

3 SOME APPROXIMATION OF THE DIELEC- 
TRIC CONSTANT TENSOR 

Now we obtain information about the electron part of the dielectric constant 
tensor e^^ (the ion part can be found in a similar manner). The polarization 
vector P is defined as 

dP ^ ^ ^ ^ 

— =j, P = Po + P'e^^[i{kf-ut)\. 

It follows that 

- ' i^^J 

where —iuP' = —eKgV, 5aj3 is the Kroneker symbol. The expression for in 

a rotating system can be derived directly from the equations of electron motion 
(the field E contains the wave field E' ~ e~*'^*) with the force (exact) 

— * — > 6 — * — * C 6 

Fe = -eE - - [vx B]+ mn'^pep + 2m[v x Q] - -B^^pep + -UpEpCz. (13) 

We suppose that the constant field B contains all components {Bop, Bq^, Bqz) 
and V ~ e"*"^*, i.e. space variations of E' are neglected. Resolving the equation 
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V = iFe/{mu) in terms of Vp,v^,Vz and extracting the ^''-dependence only (the 
constant polarization Pq and the constant part e^^^ exist), one obtains 

^ _,_mijzs-m (14) 



where = {4nK,eymy/^ = ^KJ^ + (^^j2 + j2 ^ 
(jj^e = eBQp/{mc), t^Be = e5o^/(mc) , = eBoz/{mc) - 2Q , 

here a, /? and 7 are taken from cyclic rearrangements in p, ip, z; e^^ = £*p^. The 
approach requires the following conditions: 

vt \ ka \ /lo <^ 1 , vt \ kf^ \ /<^§g <C 1 , (the smallncss of space £"-variations in 
the localization region of an electron), i.e. lo must be large enough; and, besides, 
(jj must not be near the frequency ui. It follows from the expressions (14), (15), 
that the dissipation equals zero: 

(the tensor Sap is Hermitian). The £a/3 values determine the wave polarization 

(which is elliptic in the general case) for large to. The gyration vector G can 
be found from the equation = EafjEj} + i\E x G\a , here i = \/— T. In this 
approximation the components Ga are: 



here a = p,ip, z. 

We note in the general case, that the space dependences of fields B and 
the wave amplitude need to be taken into consideration, in principle (to obtain 
rigorous expressions for £q,^ or damping of waves, for example). One may sup- 
pose, that, as a consequence of space nonhomogeneity of all characteristics, the 
rotating plasma system as a whole can be more stable with respect to exter- 
nal disturbances (the system goes from resonance influences). To impart some 
energy, the spectrum needs to be rather wide-range. 



4 THE LOCATION OF POLES 

Since the rotating plasma is a nonhomogeneous accelerated medium, the wave 
propagation has some peculiarities (damping, non-rectilinearness). In the 
general case the resonance conditions possess a relatively complex structure 
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ijj{p, (f, z,B,E,v) = 0. The derivation of the resonance conditions will be demon- 
strated below in the context of the linear kinetic description for some special 
cases. Note, that the elementary approach to the problem can be developed 
from the analysis of the exact force (13): it is necessary to find such particles 
{v{p,ip, z,t)), for which the work of the wave field A = J Fvdt is positive for 
some specific wave (in type, direction, polarization). The value A is 

A = j {-e{Ev) + mQ'^pVp + -9.p{BpVz - BzVp)}dt. 

If the mean work is positive, the Landau damping mechanism takes place. In the 
kinetic description we start from the system of self-consistent field equations: 

^ + + = 0, (16) 

rotE = — , divB = 0, 

c ot 

19^ 47r^ , ^ . 

rotB = -— — I 7 , divE = Airq, 

c ot c 

q = ej {Nifi - fe)d^p , j = ej {Nifi - fe)vd^p, 
and separate the disturbed terms (caused by the wave propagation): 

fa = f0a + 5fa, B = Bo + & , E = Eq + E' , (17) 
where a means either e or i. It is believed that 

5fa, E', B' A{p,z)expi{kf-ujt). (18) 
The undisturbed distribution function (in the rotating system) takes the form: 



K,D 



(19) 



= ^ + Ue, (20) 

where Kg is the total number of electrons in the volume V, D is the normalization 
constant, Ue is taken from (4). The expression (19) is a solution of the kinetic 
equation (16) with the force (13). Ue is the potential part of electron energy 
(related with the force F). Note that the nonconservative force ev x B is not 
appeared in the potential at all. All terms with B include the factor J7 and exist 
as a result of the system rotation only. We use the substitutions 

Eq = Eop{p, z)ep + Eoz{p, z)ez, 



9 



and k' = k + S^rjf, where 



S(rl = -i,[f'-^^} (21) 

Now we make some remarks. 

1) The accounting of the amphtude space-dependence is equivalent to including 
some imaginary part in w: lo ^ lo + iu'. For, the value co' is greater than zero 
in the case of wave damping, and it follows from (18), (21), that the appropriate 
wave amplitude A{p, z) is decreasing; the unstable case cu' < is analogous to 
wave amplitude increasing (see (18), (21)). The coUisional damping can also be 
included in the term (21). The reverse problem can be formulated: with the 
A{p,z) dependence known (using a probe wave), one should obtain information 
about collisional processes (the integral of collisions). However, it is not our 
intention to follow this way. 

2) In the general case of unsteady rotation and gravitational force existence the 
following substitutions need to be done: 

Eop Eopi = Eqp -^9p{p-, z) , sj^^J eI^^J^ = Ej^^J 'i'^^^^' 



in (20): 



where 



M{p)r 



g{p,z) = G- 



G is the gravitational constant, M{p) is the mass in the < p' < p region; in(13): 

mO^pep — > m.Q^{p, z)pep + mg{p, z). 

The derivation of the disturbed distribution function is demonstrated in Ap- 
pendix. 

Of main interest here are the poles of 6fe (they contribute to the Landau 
damping). As a consequence of different regimes of particles behaviour (see 
Appendix), the resonance conditions differ depending on the particles velocity. 
For a < 1 the term Q{(p — r) does not have any pole in an explicit form (see 
(A8) in Appendix). The imaginary part of the dielectric constant tensor arises 
from the imaginary part of S/m and from the integration of the exponent . We 
consider the case | a |^ 1, i.e. v± is large enough: 



,W ' (z) ■ 

All terms (sec (AS) in Appendix) can be expanded as a power scries in a. Inte- 
grating over the region r ^ (which is most considerable for small | a |), one 
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obtains 

Sfe = Qi'p) /o°° exp {— ■iry'ildr, where 

Yi = : 1- Aiacos((p + (po) + di + 2D - (di + D)asm((p + (fo) , 

1 + sm + cos if 

i.e. Sfe = —iQ{(p)/Yi . To find the location of poles, it is assumed, that 5 = 0; 
as a result, we have the resonance conditions: 

kzVz — uj — v±k± sm{(p + ipQ + (ps + a) 



r^{siii(^;-j + a) — 2 cos(^:-j + n) [1 + cos ^ + sin ip]} 



(22) 

(1 + cos ¥J + sin ip) ^ nfp^ + t;2 {J^ly 

Thus, the space-dependence of physical quantities leads to considerable change 
in the resonance conditions, except = (for which the well known result 
follows: kzVz — u) = 0). Note, that the resonance conditions depend on the space 
variable p and angle a. Besides, the conditions (22) are nonlinear with respect 
to v±,p>,Vz, i.e. in the general case two of these variables cannot be chosen 
arbitrarily. Therefore, the resonance region decreases. 

In spite of the existence of poles in explicit form for | a |= 1 and 62 = !> the 
latter conditions are additional; therefore, the resonances don't contribute to the 
Landau damping after integration in (A9) (see Appendix). 

Now we consider the case | a |> 1. One can easily see from the structure of 
term Q{(p — t), that there are the poles of the integrand at a > 1 for 

p — T + (po = aicsm(—-). (23) 
a 

However, for the additional Landau damping the oscillations of the exponent in 
(A8) (see Appendix) near these poles should be studied. For a > 1 the oscillations 
do not increase for (p + cpQ = — arcsin(l/a) and r = 2mr, n = 0,1,2,.... However, 
expanding Y{p,t) in r near r = 2n7r for n = 1,2,..., one can see that the 
integration in the limits from —00 to +00 in (A8) gives zero. It follows for n = 0, 
that Sfe = Q{p) exp {—iTY2}dT, where 

Ai Aiacos{p + po) {di + D)cospi 

-12 = : r : — 7 r + - 



1 + sin (/9 + COS 99 l + asm{p + po) - l{sm{p + po) + siniy^i) 

Therefore, Sfe = —iQ{p)/Y2. The third term in Y2 is negligible near p + pQ = 
— arcsin(l/a); the first term is also negligible, except p = 7r,37r/2. We suppose 
5 = and, using the definition of p\, obtain the resonance conditions: 

kzVz + k±_v±_ sin((/?3 + a + pi) - uj = Q. 

Since p is fixed, the resonance doesn't contribute to the Landau damping. The 
expression for Y2 can be used for finding approximate location of poles in the 
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case of I a |S> 1. The latter resonance gives the Landau damping. The resonance 
conditions are: 

kzVz + k±v± sm{ips + a-ip-(po)-uj = — ; — -. . (24) 

1 + cos (f + smip 

We can sec the transformation of the resonance conditions, which is a consequence 
of plasma motion and inhomogeneity. Note, that (24) gives rise to the well known 
condition kzVz — co = for k± = 0. Thus, in the general case there exist resonance 
particles (poles in Sfe), which contribute to the Landau damping. 

We make some remarks on the local-rotational description of plasma. One 
can imagine that in a variety of cases the same system is composed of different 
subsystems. The question of the scale of averaging is not trivial. Since plasma 
characteristics depend on the regime in it, it is clear that the averaging over 
different local subsystems can lead to different results. Here the question arises, 
what particular averaging is correct? It is possible that this case just represents 
the local-rotational regime (since there is additional damping as waves propa- 
gate in an accelerated medium). The basic ideas of the description consist in 
the following. Instead of using averaged plasma characteristics for finding some 
physical quantity, we shall seek the quantity in a small rotating subsystem with 
definite parameters (any plasma particles are involved in the local rotational 
movements). The net result (the specific physical quantity under investigation 
in a plasma system as a whole) can be obtained by averaging the quantity over 
all possible directions of subsystem rotation. Actually, we can imagine that the 
plasma (even though it is collisionless) consists of rotating parts with radius 
r ^ vd and linear velocity v ^ Vm, where r^) is the Debye radius 



l.^^^^^K„Mi . (25) 

Vm is the most probable velocity 

Vm = \— . (26) 
V m 

Therefore, the characteristic frequency of the local rotation is 

^l-^ (27) 

(collisions produce local rotations with r < y , where /q is the free path length). 

In the general case the solution (but not equations) for a physical quantity 
under study (transport coefficients, for example) should be averaged over the 
sequence of local rotating "cells". For a "cell" the kinetic equation takes the 
form (for electrons): 

df df ^ e - 

-^+v^ + {-eE - -[v X B] + 2m[v X ni]- 
ot or c 
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^[B,e,-B,e,]+mi}fp)^ = ^^^ . (28) 

C op T 

This way of inclusion of local regimes (with subsequent averaging over all 
directions of Cli) can be adaptable to the strong turbulence. 



5 THE SELF-SUSTAINED ROTATION 

First of all, we note that any rotation induces some additional electric field Eq, 
which is directed towards the axis of rotation (z-axis) . Assume that the initially 
non-rotating plasma system resides in the external magnetic field Bq = e^iJo- 

1) If the electron surplus exists (A > 0), then the radial field Eop is directed to 
the 2;-axis. The Eop x B drift of charged particles results in the movement round 
a circle. The originated centrifugal force causes partial separation of charges. 
Therefore, the field | Eqp \ will increase resulting in sustaining of rotation. Thus, 
a collective rotation with O directed as B (parallel) arises spontaneously in a 
plasma system with an electron surplus. 

2) If the ion surplus exists (i.e. A < 0), then the radial field Eqp is directed 
from the z axis. The E^p x B drift initiates rotation with being antiparallel to 
B. However, the partial separation of charges, which is caused by a centrifugal 
force, results in decreasing | E^p \. This effect counteracts the rotation. Since 
there exists the rotation with O antiparallel to S in a positive (only) particle 
system (it is some limit), there exists some critical value of the ion surplus Ac 
(A = Ac < 0), starting from which the spontaneous rotation of plasma (fi an- 
tiparallel to B) is possible. 

We shall show that the spontaneous rotation with Q parallel to B can arise in 
an initially neutral plasma system (A = 0). The condition, under which the self- 
rotation arises, is QRq < Vd\ here Vd is the drift velocity, i.e. BqVLRq/c < —Eop. 
The equality in the condition is possible for J7 = and Q, = ^self- Which equi- 
librium is stable? We note, that the equiprobable rotational small disturbances 
(fi 0; or, factually, the ratio of the rotational energy to kT is small) with 
Q parallel and antiparallel to B are not equivalent for the plasma system as a 
whole. As a result, plasma will rotate on the average in some definite direction: 
any rotational disturbance produces the radial field Eop directed toward the ax:is 
of rotation, but the appeared Eop x B drift promotes rotational disturbances with 
Q parallel to B only. We have for the other limit (the case of externally initiated 
large Q,, i.e. MjQ^p^/(2fer) 3> 1) the inequality Opo > Vd, since the limitations 
on Eop always exist for any plasma system: — oo < £^op'" — ^Op < Eg^"^ < oo. 
Therefore, the drift opposes the rotation. Thus, the rotation tends to some equi- 
librium ^Iseif 7^ for A = 0. 

In the general case, when the rotation (ri is parallel to Bq) arises by any 
reason, the separation of charges and electric field arise. As a result, the E x B 
drift of charged particles (the movement round a circle) arises. This movement 
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will either add to the initial rotation or subtract from it. Eventually, the stable 
rotation is established. Prom the equation 

Vd = ^{p, z)p, 

where the drift velocity is 

_ _ cEop{p, z)e^ 

Vd — r> 
-DO 

it follows, that the stable rotation frequency is 

^self{p,z) = -'-^^, (29) 
Pm 

with Eop{p,z) taken from (6) and | |^| I- 

Under an assumption, that the originated magnetic field Bz is comparable to 
Bq, the following substitution needs to be done in equation (25): Bq —> Bq + B^ 
with Bz taken from (10). In this case the self-sustained rotation of plasma takes 
place. 

We make some remarks. 

1) Indeed, to evaluate the rotational state stability one can involve the energy 
balance. In this case the energy surplus characterizes qualitatively the rate (or 
time) of transition. However, a rather complex integral system of equations for 
^seif can be solved numerically (or by a perturbative method) for all specific 
values of plasma system parameters only. 

2) ^self tends to zero for the temperature high enough. 

3) If the external magnetic field Bq is zero, then the self-sustained rotation is 
impossible (the velocity p^lgel f must be greater than the speed of light) . 

4) If several quasistable field distributions (distinguished by a potential barrier) 
can exist, then the existence of several quasistable rotating regimes (with ^l^if) 
is possible. 

6 CONCLUSIONS 

The generation of electric and magnetic nonhomogeneous self-fields in rotating 
plasma can be explained by various actions of the centrifugal force on dissimilar 
particles. There are different /o-distributions of these particles and the partial sep- 
aration of charges. Therefore, the electric field arises (sec (6), (8)). The magnetic 
field (see (9), (10)) exists as a result of rotation of charged regions (the mag- 
netic actions are not compensated). As for relation to the physics topics from 
Sec. I, there are new peculiarities regarding the self-consistent description. All 
self-fields Ep{p, z), Ez{p, z), Bp{p, z), Bz{p, z) are taken into consideration. The 
self-sustained rotation (which origins from initially neutral state) exists as a re- 
sult of the particles drifts in the electric self-field and external magnetic field (see 
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(29)). This is a new result. The dielectric constant tensor in rotating plasma 
differs from that in non-rotating plasma. As the result of plasma inhomogeneity 
and accelerated motion, the resonance conditions change (see (22), (24)). There 
are new poles in a disturbed distribution function and, therefore, new terms in 
the Landau damping. The local-rotational description (over rotating "cells") of 
plasma can be useful. 

APPENDIX: Derivation of the disturbed 
distribution function 

The problem is to obtain the disturbed electron distribution function. We 
substitute (17-21) and the force (13) into the kinetic equation (16), correct to 
the second order {df,E',B' are small). Then it turns out that 



where 



mn^ip,z)p-^ --[vx Boi]-^ - eE'J-^ = 0, (^1) 



mp mc e 



Eq^ — Eqz BopCz 



c 

Transferring to the cylindrical coordinate system (vz, v±, cp) in the velocity space 
(Fig. 2), one gets 

i[kzVz + k±v± cos — Lo + S{'r)pv± cos{(p — a) + S{r)zvz]5fe H — -^{vE')— 
-{B Vz - B^v± cos [tp - a)) + |i2ipcos(<^ - a) - sm {(p - a)vz\- 



[^Be - — sm((^-a) -u;^gCos(9?-a)— J- 



v±_ v±_ dip 



sin - aW - = , (^2) 



d 

to the ray. Then we obtain the equation 



where is the given value, sin a = ^, d is the minimum distance from the center 



d5fe aisimp + bicosip + di 

+ ^ ■■ —r 7-r°J^ = '5(W, (^3) 

oip a2 sm ip + 02 cos ^p + \ 

where 

S{f)v_i_d k±v± + S{f)^J fP' - (Pv± kzVz - OJ + S{r)VzZ 

'^i = — W~ ' " ¥) ' " ¥) ' 

^Be ^Be ^Be 
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FIG. 2. The cylindrical coordinate system in the velocity space. 



z 



nip . J^lv^ njp J^lv, . 

02 = TT Sin a Py — cos a , 02 = rr cos a pi — sm a 

(2) (2) ' ^ (2) (2) 

<3(w = -fiy- -> ^ = — ^{(^^^) i^p^- 

a;^g(a2sin(/? + 62COS(^+ 1) ^ 

cos - a)]} + - sin - a^^- 

[VlIp cos (99 - a) - Wsgi;^ sin - «)] |^ • 
The solution of equation (A3) is 

Jc 



where 



that is, 



dy ai sin (p + bi cos if + di 
dip a2 sin + 62 cos ip + 1 ' 

02 sin + 62 cos (p + 1 



bia2 - aib2 , 



1 + tan ^ 



+ 



6162 + 0102 6162 + 0102 \ P dp' 

02 + ^2 a| + 6| ^ y a2 sin y?' + 62 cos y?' + 1 ' 

dip' 



a2 sin + 62 cos (p' + 1 
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: arctan 



(1 - 62) tan 1+02 



O2 0.2 



bl - a\ 



2) 



:ln 



(1 - 62) tan ^ + 02 



62 + 



^Jhl + al - 1 (1 - 62) tan f + 02 + + ai - 1 



, 6i + ai > 1, 



3) 



In 



02 



1 + 02 tan 



4) 



«i + 6i 



1. 



02 + (1 - 62) tan I 

Under the (/9-periodicity condition on 5fe it turns out from (A4) (with the substi- 
tution = (/? — r), that this expression (note that all functions are (/^-periodic) 
can possess the 99-periodicity for C = ±00; the value C = 00 obeys all transitions 
(in limits) to the well-known results only; therefore, 



5fe 



o*[y(v-T)-2/(¥')] 



Q{}P - T)dT. 



The substitutions 
62 



(/5o = arctan ■ 



vi_ [ki_ sin (993 + a) + pS sin (^3] 



02 



aw 



(^) 

Be 



tan (^3 



+ 



(p)n2,,2 



Bey 



D = 



v± [k± cos {a + (fs) + pS cos (^3] 



Ci = di + -D , = arctan (l/Va^ — 1) 
and some rearrangements give: 



aa; 



(^) 

Be 



y(^,r) =Ai hi 
where Wi = 

1) 



;>oo 

Sfe= ex.p{iY{ip,T)}Q{ip -T)dT 
Jo 

[1 + a sin ((/5 + (^0 - t)] [1 + tan ^] 



{A8) 



[l + asm{ip + ipo)] [i+tan(^)] 



; arctan 



sm ■ 



2) 



1 



3) 



:ln 



sm 



cos 5 + asin((/p + (^0 — §) 



-Dr + CiVFi 



o |< 1 



sin(£±i^±£i)cos( 



2 ' 



02 



1 + 02 tan ^^2^ 



1 + 02 tan ^ 
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2 sin I 



cos I + sm{(p + (^0 - i) 



1 + asin((^ — T + (^o) 
— S^uj, cos (y? — a — r))] + av± sin((^3 + a — + r 



d5f, 



dv± 



/e^^i (P) ■ / x\ 
■ (^-^ - <^BeV± sm(<^ - a - t)J 



where the vectors ^" , B" are the vectors E', B' turned at angle r in the plane 
normal to Bq. 

Note that expression (A8) is rigorous. As the first approximation, one can 

use 

d8f{ip - r) _ dSfuj^ - r) d5f{ip - r) _ d5fMi<f - r) 
dvj_ dv± ' dvz dvz ' 

where the disturbed part of the Maxwell velocity distribution function is 

roo ^1 

Sfui^) = / exp {-iair' - 2iP cos{ip - — ) sin —)Q*{^ - T')dT', 
Jo 2 2 

kzVz-(^ r, k±v± e{vE')fM 

^Be ^Be ^ ^Be 

Here is the Maxwell velocity distribution function. Note from (AS), that 
there exist different regimes of particle behaviour (they depend on the velocity 
of particles). The dielectric constant tensor can be found from 



-Ea = — / / VaSfedVzd(pV_idv± . {A9) 

lUJ Jo Jo J-Qo 



'^ap '^ap J-, c 

Ep = — I I 

47r lio Jo Jo 

The dielectric constant tensor differs from that in non-rotating plasma. 
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